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(54) Method of correcting a radiation image for defects in the recording member 



(57) Defects in a member for recording a radiation 
image are corrected by generating a calibration matrix 
comprising for each pixel a correction value and by 
applying corresponding elements of the calibration 
matrix to the pixels of a recorded image. The values of 
the calibration matrix represent in each pixel the devia- 
tion of the value that is effectively obtained in a pixel 



when the recording member is subjected to a flat field 
exposure and the value that would be expected in the 
pixel. Different methods are described for determining 
the expected value in a pixel. 

Further, an overlay image is generated for a recording 
member wherein local defects are visualised. 
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Description 

1. Field of the invention. 

5 The present invention is in the field of medical imaging. The invention more specifically relates to the detection and 

locating of defects of any type and size in a recording member for storing a radiation image. The invention further 
relates to the correction of such defects in a radiation image. 

2. Description of the state of the art. 

TO 

In the field of digital medical imaging a wide variety of image acquisition techniques have been developed that 
render a digital representation of a medical image. Among such techniques are computerised tomography, nuclear 
magnetic resonance, ultrasound, detection of a radiation image by means of a CCD sensor or a video camera, radio- 
graphic film scanning, techniques using Selenium and a read-out system such as thin film transistors etc. 

15 In still another technique, X-rays transmitted by an exposed object (such as the body of a patient) are stored in a 
photostimulable phosphor screen. Such a photostimulable phosphor screen comprises a layer of photostimulable iumi- 
nescent material which comprises a phosphor, for example a europium-activated barium f luoro-halide, and a binder. 
The phosphor has the characteristic that it can be energised to an excited state by X-rays, and can then be stimulated 
by light within a first wavelength range to return to the ground state with the emission of light within a second wavelength 

20 range. The stimulating radiation is arranged to have a different wavelength from the emitted light. With regard to the 
specific example wherein a radiation image is stored in a screen comprising a europium-activated barium f luoro-halide 
phosphor, the stimulating light is situated within the range of 600-700 nm and the emitted light is situated within the 
range of 350-450 nm. 

In order to read the image stored in an exposed screen, the screen is transported past a scanning station where it 
25 is two-dimenstionnally scanned with stimulating light. For this purpose a light beam of a wavelength within the stimulat- 
ing wavelength range of the phosphor is deflected so as to scan the screen in a first direction (called main scan direc- 
tion) while it is transported in a second direction that is perpendicular to the main scan direction. 

Light emitted upon stimulation is guided through a light guide member onto an opto-electric transducer such as a 
photomultiplier that converts the emitted light into a corresponding analog electric signal. This analog electric signal is 
30 then converted into a digital image signal by means of an A/D converter. 

A digital image representation acquired through the use of a photostimulable phosphor screen as an intermediary 
storage means or a digital image representation acquired by any other kind of detector and/or image acquisition sys- 
tem, always comprises noise to some extent. 

Especially in the field of medical imaging wherein systems are developed for extracting information from a patient's 
35 body and for presenting this information to a physician to allow him to perform a diagnosis, it is extremely important that 
all imaging defects in the image are detected and signalled since they may interfere with relevant signal components 
and may as such have an influence on the diagnosis. 

The present invention deals in particular with the detection of defects on the digital image representation of a med- 
ical radiation image. Still more specifically the present invention deals with defects on a radiation image that is stored 
40 in a photostimulable phosphor screen that is read out according to the above described read-out process. 

Defects on a stimulable phosphor may originate from several causes. In the context of the present invention, these 
defects may be classified into three classes: (1) point defects originating from grain size impurities in the moulding proc- 
ess of the phosphor screen, or from dust during subsequent use. The spatial dimension of point defects is smali with 
respect to both principal image dimensions, their gray value amplitude exceeds the imaging chain noise amplitude in a 
45 case-dependent extent. (2) line-like or curve-like defects, originating either from smear during the production process 
or from mechanical wear during use. These class of defects may further be subdivided by (a) pure line defects such as 
- but not limited to - produced by mechanical rollers; (b) scratches due to e.g. handling and damage; (c) curve-like 
defects due to any kind of one-dimensional distortion process. All of these subclasses of linear defects has large ratio 
of length along the axis of the defect w.r.t. the across dimension; the amplitude can vary widely, the weakest amplitude 
so (as produced by e.g. smear) being considered unacceptable. (3) area defects such as spots, and global impurities e.g. 
coming from the moulding process at production or moisture during use. This class of defects shows dimensions in 
either of the principal image axes. Still further causes of sources of defects may be thought of, but the shape of any of 
the resulting defects can always be classified as either a point, a line or an area defect. A common characteristic shared 
by of all of the aforementioned classes of defects is that they may be present already at fabrication time or come into 
55 being during use of the imaging plate. This fact necessitates the periodical inspection of imaging plates to detect the 
presence of (newly) appeared defects, at regular intervals of the life time of the plate. 

The early interception of these defects and the subsequent corrective action is indicated at least in the production 
stage and is particularly motivated during use in critical application areas such as mammography. Regular screening of 
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the total installed base of imaging plates as to the presence of any of the mentioned defects is mandatory. This kind of 
quality control is commonly known as visual inspection. Human inspectors, however, are subject to subjective percep- 
tion and assessment. Moreover, the presence of subtle abnormalities such as dust, small scratches and global impuri- 
ties may be poorly detectable or even remain unnoticed. However, these defects leave a clearly visible imprint in the 

s digital image, and may severely impair correct diagnosis by the radiologist. Dust and scratches could possibly be con- 
fused with potential microcalcif ications on digitized or digital mammograms. Mammography is a typical application area 
requiring the highest image quality standards available. The field of research aiming at developing automatic detection 
methods for any of these defects by computerized inspection is commonly known as automatic visual inspection. 
The result of the automatic visual inspection may either be (1) the complete rejection of an imaging plate due to an 

10 unacceptable defect, according to a quantifiable measure or (2) the marginal acceptance of the imaging plate with sig- 
nalling and indication of zones requiring special attention or (3) acceptance of the image plate. 

Automatic visual inspection is a well developed field in digital image processing, such as presented in the F. van der 
Heijden, Image based measurement systems, Wiley, 1994. However, very little mention is made in the literature of 
methods applicable to the highly specialized area of film inspection and phosphor imaging plate inspection, in P. 

15 Dewaele et a/., An application of the Lily software package to defect inspection in unexposed radiographic film, Proc. of 
the SPIE, The Sixth Meeting in Israel on Optical Engineering, Vol. 1038, Dec. 1988, pure one-dimensional techniques 
based on edge and level measures are developed for fast on-line inspection of radiographic film. However, they assume 
linear defects in the along direction only. Automated inspection of phosphor plate has to cope with a much wider variety 
of defect geometry, necessitating the use of two-dimensional techniques.Currently, inspection of phosphor imaging 

20 plates is done manually and on a representative sample of a production batch. However, in view of the benefits given 
before, an automatic inspection is highly desirable. 

3. Objects of the present invention. 

25 It is a primary object of the present invention to provide a method to detect and locate point, line and area defects 

of a recording member in which a radiation image is stored, and to assess the degree of severity, to delineate the region 
of support of the defect and to correct an actual radiation image for these defects. 

It is a further object to provide an accompanying overlay image showing the defective zones on the recording mem- 
ber requiring specific clinicians attention, and to provide a two-dimensional calibration matrix supplying correction fac- 
30 tors to correct a registered actual image on a pixel -by -pixel basis. 

Still further objects will become clear from the description given hereafter. 

4. Summary of the invention. 

35 The objects of the present invention are achieved by a method of correcting a radiation image stored in a recording 
member for recording member defects comprising: 

generating a calibration matrix for a recording member, said calibration matrix comprising for each pixel on said 
recording member a correction value that represents the deviation of the value that is effectively obtained in a pixel 
40 and a value that would be expected, when said recording medium is subjected to a flat field exposure and subse- 
quently read-out by means of a calibrated read-out system, 

generating a digital signal representation of the radiation image by reading the image by means of a calibrated 
read-out system, 

correcting each pixel of said radiation image by applying to a pixel a corresponding value of said calibration matrix. 

45 

A calibration matrix is generated as follows : 

first the recording member is exposed to a flat field exposure of radiation, 

next the image stored in the recording member that has been exposed to the flat field exposure is read-out with a 
so calibrated read-out system (so that influences on the read-out signal that do not originate from a defect in the read- 
out member are eliminated) 

next pixel values of the flat field image are averaged to generate an average value of p R , 
then a smooth background image 9r(U) *s calculated by one of the methods described below, 
finally the values of the calibration matrix are calculated either as UR-g R (i,j) in the additive case wherein pixel values 
55 of the radiation image are corrected by adding to a pixel value a corresponding value of the calibration matrix, or 
as *iR/gR(i,j) in the multiplicative case wherein pixel values of the radiation image are corrected by multiplying a 
pixel value of the radiation image with a corresponding value of the calibration matrix. 
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The generation of the calibration matrix and in particular alternative methods for the generation of the background 
image will be explained below in more detail. 

In the context of the present invention the term 'radiation' refers to penetrating radiation of any kind such as X-rays, 
gamma rays etc. 

5 A recording member in the context of the present invention is a member wherein a radiation image, for example an 

X-ray image, can be temporarily or definitely stored. Different types of recording media may be used. One example of 
a reusable recording member is a photostimulable phosphor screen. An example of a non-reusable member is a radi- 
ographic film. Still others may be envisaged, such as different screen-film combinations. 

By means of the term 'flat field exposure' is meant a uniform exposure of the recording member to penetrating radi- 
10 ation so that ideally in every pixel of a defectless recording member an identical value would be read-out when using 
an optimally calibrated read-out system. 

A calibrated read-out system is a system that has been corrected for any non-uniformities in the read-out signal that 
originate from the system's components itself. Causes for such non-uniformities can for example be attributed to the 
means for directing light onto an opto-electric transducer. These phenomena as well as examples of adequate correc- 
ts tion methods have been described in extenso in the state of the art literature, for example in US patent 4,885,467 and 
in Research Disclosure N 352, page 484-485 (1993). 

The term 'corresponding' refers to corresponding elements ifTthe actuarimage and the calibration matrixrboth- 
being brought in the register as explained furtheron. 

20 5. Brief description of the drawings 

Particular embodiments of the invention will now be explained in more detail with reference to the accompanying 
drawings wherein 

25 Figure 1 generally illustrates an apparatus for reading out a photostimulable phosphor screen, 
Figure 2 schematically illustrates the method of the present invention, 

Figure 3a depicts a typical scan line of a stimulable phosphor screen subjected to a flat field exposure, 
Figure 3b shows the same scan line having been subjected to background removal, 
Figure 3c shows the result of a thresholding operation performed on the scan line of Figure 3b, 
30 Figure 4 depicts a particular embodiment of a hierarchical multi-resolution image decomposition method and the 
generation of a background image starting from the decomposed image, 
Figure 5 shows an example of a filter used in said decomposition method, 
Figure 6 illustrates the calibration matrix computation and the error detection, 
Figure 7 illustrates the correction of an actual image. 

35 

The method of the present invention will be explained in general with reference to figure 2. The part relating to the 
acquisition of the digital signal representation of a radiation image is illustrated in figure 1 for a specific example wherein 
a radiation image is stored in a photostimulable phosphor screen. 

After a stimulable phosphor sheet 1 is exposed to radiation, such as to X-rays, passing through an object to have 
40 a radiation image stored thereon, it is sent to the read-out station shown in Figure 1. A laser beam 3, having a wave- 
length of 633 nm. is emitted from a helium-neon laser source 2, and directed towards a galvanometer mirror 8. Drive 
means 6 cause the galvanometer mirror to reciprocate in a triangular wave pattern. A light chopper 4, including a rotat- 
ing disc segment 5, is positioned to enable the disc segment to interrupt the laser beam pathway during the galvanom- 
eter retrace step. 

45 Various laser beam focusing devices, known in the art, not shown in the drawings, ensure a uniform beam diameter 
during scanning of the beam on the phosphor sheet and also ensure that the uniform angular velocity of the reciprocat- 
ing mirror 8 results in the laser spot travelling across the phosphor sheet at a uniform linear speed. The laser beam 3 
is one-dimensionally deflected by the galvanometer mirror 8 and by a plane reflection mirror 9 onto the phosphor sheet 
1 . Transport means, not shown, are provided to transport the sheet in the direction of the arrow 1 1 , to enable the whole 

so sheet to be scanned in a uniform manner. 

Positioned closed to, but behind the scanning line of the laser beam 3 on the phosphor sheet 1 , is a light guide 12 
which receives light emitted from the phosphor sheet 1 but is shielded from direct exposure to the laser beam 3. The 
light guide 12 consists of individual optical fibers that are bundled at the output end. The output end of the light guide is 
positioned adjacent a photo-detector 13, which produces an electrical signal dependent upon the light intensity falling 

55 thereon. 

Suitable electric connections are made to pass the output signal from the photo-detector 13 to a computer 20 which 
serves to control the light chopper 4 and the galvanometer mirror drive 6 and to perform the method of the present 
invention as well as other kinds of image processing techniques on the output signal. 
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The output of the computer is additionally connected to a display unit 21 , such as a VDU screen. Alternatively or addi- 
tionally the output of the computer 20 is used to produce a permanent record of the image. The output can also be fed 
to a workstation for off-line processing and/or to an archive station etc. (not shown). 

The read-out signal obtained at the output of the photodetector 13 was fed to a square root amplifier 30 where it 
5 was amplified. It was then sampled by a sample and hold circuit 31 , and converted into a 12 bit signal by means of an 
analog to digital convertor 32. 

A shading correction 40 was then applied to the read-out signal by multiplying each pixel value with an appropriate 
correction value determined in advance. 

The digital signals thus obtained were sent via switch 33 either to a section 34 where correction values were cal- 
io culated according to the method of the present invention, or to a noise reducing correction section 36 where a read-out 
signal was corrected using these correction values. 

The procedure for calculating correction values is as follows. 

Before read-out of actual diagnostic images was performed, each screen that was to be used in the read-out appa- 
ratus was subjected to a uniform flat-field exposure (in the absence of any patient or X-ray modulating object). 
is The uniformly exposed screens were then consecutively fed into the read-out apparatus for identification and read- 
out. 

Identification was performed by reading the identification information stored in an electronic memory (an EEPROM) 
provided on the cassette conveying the photostimulable phosphor screen. 
Image read-out was performed following the method described higher. 
20 The read-out signals were then subjected to square root amplification, sampled and converted into a digital signal 
as has been described higher. 

The digital signals were then applied via switch 33 to the correction value calculation circuit 34 where for each individual 
pixel of a screen a corresponding correction value was calculated as will be explained extensively below. 

Then, the matrix of correction values for each pixel in a screen, in the following called a calibration matrix, was 
25 stored in memory 35 together with the associated identification item. Typically about 30 different photostimulable phos- 
phor screens are used in a single read-out apparatus. 

This procedure was performed for every photostimulable phosphor screen that was to be used in the read-out appara- 
tus. 

Instead of storing calibration matrices, one could perform the procedure of determining a calibration matrix before 
30 every exposure of a recording member to a radiation image. Although being extremely accurate, this procedure would 
be highly time-consuming. 

It is also possible to perform the procedure according to the present invention only for these photostimulable phos- 
phor screens that are used for specific examination types such as e.g. mammography in which confusing of e.g. micro- 
calcification with a defect on a recording member could easily occur. This would have detrimendous effects on the 
35 diagnosis. 

During normal use of the screens, a radiation image (of a patient or of an object) stored in a specific screen was 
identified and read-out and the digital image signals representing the radiation image were fed via switch 33 to correc- 
tion circuit 36. The noise reduction circuit then retrieves from memory 35 the correction values pertaining to the identi- 
fied screen and carries out correction by addition or multiplication of read-out pixel values with/or the corresponding 
40 correction values. 

Next, the corrected signal was applied to a processing circuit 37 where it was further subjected to various kinds of 
processing such as frequency and gradation processing etc. which methods fall out of the context of the present appli- 
cation. 

Finally a hard or soft copy image was generated on an output device 38 (a hard copy recorder, display unit etc.). 
45 As will be explained below, local defects that are detected after background removal are threshold with a floating 
threshold value T in circuit 34 (see Figures 3b and 3c). The result of this thresholding operation is a digital overlay image 
(bit map image) which is stored in memory 35 and retrieved at the time an actual image stored in an identified photo- 
stimulable phosphor screen is corrected. 

The digital overlay image is applied to a hard copy recorder 38 for the production of an overlay image which represents 
so defective area as coloured regions or contours. 

The overlay image is superimposed by the radiologist on the hard copy or soft copy of the actual radiographic image so 
that he is informed of the location of such local defects and will not confuse these defects with elements that are relevant 
with regard to the diagnosis. 

The digital overlay image is further applied to a characterisation and classification circuit 39 wherein the defects are 
55 classified according to a number of characteristics as explained in detail furtheron. 

In this way a quality report of the recording member can be created and output. On the basis of this report acceptance 
or rejectance of a recording member can be decided. 

The steps of the method of the present invention that have been described in general will be explained below in 
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more detail. First a number of procedures are explained for calculating the background image g R (i.j) used for generating 
a calibration matrix. 

4. 1. Modeling of the defectiess recording member 

5 

When the recording member is exposed to a flat field exposure, the signal is ideally expected to be a constant with 
noise superimposed. The noise consists of Poisson distributed quantum noise, dependent on X-ray exposure and fixed 
structured noise, independent of X-ray exposure. At moderate amount of exposure, the first noise source dominates, 
and by subsequent square root compression, it can be modelled as independent additive Gaussian noise with predict- 
10 able standard deviation. The signal is thus ideally composed of a spatially constant mean value with signal independent 
additive noise. 

As enumerated in the foregoing, the ideal model is distorted according to a variety of defect sources. Fig. 3a depicts 
a typical scan line of a stimulated phosphor subjected to a flat field exposure. In contrast to being a constant, it shows 
a slowly varying background immersed in noise. The defectiess member is characterized by deriving an analytical or 
15 discrete function representing the background latitude. 

Several methods have been tested to model the background. They are either analytical or empirical in nature but 
all have in common that they are level-based methods andiwill entail I evel -based def ec^ dete^ionrThe sequershortty 
describes these methods. 

20 4.1.1. Global 1- and 2-dimensional regression 

Regression is a widely employed technique and ready-to-use implementations are available such as presented in 
Press. W.H. et a/., Numerical Recipes in C, 2nd Edition, Cambridge University Press, 1992, Chapter 15: Modeling of 
Data. Regression is a method where the N=K.L data points, in casu the 2-dimensional array of gray values, are fitted 
25 to a model that has M adjustable parameters a,, j '= 1 .. M . By the qualifier 'global' is meant that the spatial extent of the 
data covers the complete image, that is K and L denote the image dimensions. The error minimized is the least-squares 
error 

/=1/=1 



where the z } j denote the actual gray values, z(/, j\ a v ..a M ) denote the functional relationship between independent and 
35 dependent variables and the ct/j are the standard deviation of the measurement error on z t j. Least squares fitting is a 
maximum likelihood estimation of the fitted parameters when the a/j are constant.This maximum likelikhood estimate 
amounts to a chi-square fitting when cjjj are non-constant and supplied for each gray value. By taking the derivative with 
respect to the parameters a k , we obtain equations that must hold at the chi-square minimum, 

o-iz(^)(^fer^ -i « 



45 This is in general a set of M nonlinear equations for the M unknown a k . 

The solution of the system of equations becomes mathematically more tractable when the model is a linear com- 
bination of any M specified functions of the spatial coordinates /, j. This kind of regression is termed as general linear 
least squares and has the following form of the model 



50 



M 



55 where the basis functions Z 1 (/, j) Z M (i, j) are arbitrary fixed functions (but possibly non-linear functions of the spatial 

coordinates x, y). The computation has further been speeded up considerably by randomly sampling a fraction of the 
total number of data points in the image plane either by using random number generators or by quasi-random number 
generators such as Sobol* sequences, introduced in e.g. H. Press et a/., Numerical Recipes in C, Cambridge University 
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Press. 1992, Chapter 7: Random Numbers. 

A special case has been tested using polynomials as the basis functions, it is easily conceivable that any other suit- 
able basis functions could be used without departing from the general principle of least squares fitting. By way of exam- 
ple, but not restricted to it, the polynomial up to exponential degree 4 has the following form 

5 

y) = 3^ + a 2 x + a 3 y + a 4 x 2 + a 5 y 2 + a 6 xy + a 7 x 2 y + a 8 xy 2 + a 9 x 2 y 2 

representing a case with M = 9. 

In a further refinement to speed up computation without compromising accuracy, the model of the background gray 
io value variation has been simplified to the pure one-dimensional case, whereby each scan line is modelled as a polyno- 
mial of degree M - 1 : 

2 M 1 

z(x, y) = a t + a 2 x + a 3 x +...+a M x 
15 Again here, the basis functions can be any linear or non-linear function of the place domain coordinates. 

. 4. 1.2. Local 1- and 2-dimensional spline fitting 

The picture can be thought of as a surface by using the brightness value as the third coordinate. The surface can 
20 be too complex however to be described by a single equation over the domain of the image, such as computed by global 
regression, and therefore the need for piecewise approximation arises. Such piecewise representation is offered by 
splines such as presented in T. Pavlidis, Algorithms for graphics and image processing, Computer Science Press, 1982, 
Chapter 1 1 Curve fitting with splines, Chapter 1 2 Approximation of curves and Chapter 1 3 Surface fitting. According to 
the nature of representation, splines can be subdivided into either interpolating or approximating, having either fixed 
25 knots or variable knots. Subdividing the image into tiles of equal size (200x200 pixels), the fitting problem belongs to 
the class of approximating splines with fixed knots. Knots (or control points) are held fixed here, since increasing the 
number of knots in areas of higher local variance is undesirable, precisely because they may represent defects. 
In contrast to the discrete methods used hereafter, the analytical representation enables much more compact storage 
of the calibration matrix, to be defined next. Instead of storing the individual additive or multiplicative correction coeffi- 
30 cient for each pixel, the coefficients of the approximating splines are stored. The actual background values are recon- 
structed at the time of processing the image. Although the storage capacity needed and the computation time required 
for spline fitting are larger than those of the global regression method, the benefit offered by splines is that they can 
more accurately model the background variation. The computation has been speeded up in an alternative implementa- 
tion by computing the approximating splines for successive one-dimensional scan line profiles. 

35 

4. 1.3.Background determination by smoothing 

The former methods are analytical, that is, they express the background as a linear combination of analytical func- 
tions. A discrete representation offers the advantage that no prior knowledge as to the nature of the background varia- 

40 tion is required, so they can handle much more complex variation. On the other hand, the storage capacity needed for 
each calibration matrix is proportional to the number of pixels. A local discrete approximation of the gray value surface 
can be obtained by smoothing with a smoothing kernel of appropriate size. The values of the coefficients of the smooth- 
ing kernel are well documented in the open literature and can typically represent a box filter with constant coefficients 
or a filter with coefficients derived from a gaussian: W.K. Pratt, Digital Image Processing, 2nd ed., Wiley Interscience. 

45 These kernels are separable and can be obtained by successive horizontal and vertical convolution. The box filter has 
the additional property that it was implemented in a recursive way, using the previous smoothed value, and one addition 
and one subtraction per direction and per pixel. Smoothing filters are convolutions and can also be obtained by multi- 
plication with the kernel's transfer function in the Fourier domain. 

so Formulas implemented 

2 dimensional convolution with an odd-size {2M + 1)x(2A/ + 1) kernel 

M N 

55 ~z{x.y)= ]T £ h(m,n)z(x-m. y - n) 

m=-Mn=-N 
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with 



10 



h(m J n) = 



1 

(2M + 1)(2N + 1) 



1 1 

1 ... 
1 



1 1 1 



1 ... 1 



15 



20 



25 



for the box filter and 



2 2 M N 

h{ mrn)^r-A-exp{ ^^with A s u ch th at ^ ^-/7<m-/7)-=-1- 
2<j 



m=-Mn=-N 



for the gaussianiy shaped filter. Typical size used is M = N = 65 . 
Separable implementation on rows and columns 



M 



r(x,y)= ^ h r {m)z(x - m, y) 



30 



~z(x, y) = c(x, y) = X h d n ) r ( x > Y ~ n ) 
n=-N 



with 



35 



40 



45 



50 



h w = fab) 11 1 -• 1] h ° (m) = (drr) 11 1 - 11 

for the box filter and 

h r (m) = 4, «p(^i) M"> = A c exp (^) 
for the gaussianiy shaped filter with A r and A c such that 

M N 

£ 21 Mm.")-!. For M = N,A r = A c = VA. 



m=Mn=-N 



Recursive and separable box smoothing implementation after initialization 
55 r(x, y) = r(x - 1, y) + z(x + M, y) - z(x - /W - 1. y) 

i(x, y) = c(x, y) = c(x, y - 1) + r(x, y + A/) - r(x, y - A/ - 1) 
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4. 1.4. Background determination by hierarchical decomposition and reconstruction 

A common framework for both local and global distortion modeling on a discrete grid, is obtained by a pyramidal 
decomposition and reconstruction such as disclosed in EP 0 527 525. A preferred embodiment of the decomposition 

5 process is depicted in fig. 4. The original image z(/, j) is filtered by means of a low pass filter 41 , and subsampled by a 
factor of two, which is implemented by computing the resulting low resolution approximation image only at every 
other pixel position of every alternate row. A detail image r 0 at the finest level is obtained by interpolating the low reso- 
lution approximation with doubling of the number of rows and columns, and pixelwise subtracting the interpolated 
image from the original image z. The interpolation is effectuated by the interpolator 42, which inserts a column of zero 

w values every other column, and a row of zero values every other row respectively, and next convolves the extended 
image with a low pass filter. The subtraction is done by the adder 43. The same process is repeated on the low resolu- 
tion approximation instead of the original image z, yielding an approximation of still lower resolution g 2 and a detail 
image r-, . A sequence of detail images r h / = 0..L - 1 and a low resolution approximation g L are obtained by iterating the 
above process L times. The finest detail image r 0 has the same size as the original image. The next coarser detail 

15 image ^ has only half as many rows and columns as the first detail image r 0 . At each step of the iteration the maximal 
spatial frequency of the resulting detail image is only half that of the previous finer detail image, and also the number of 
rows and columns is halved, in accordance with the Nyquist criterion. After the last iteration an image g L is left which 
can be considered to be a very low resolution approximation of the original image. In the extreme case it consists of 
only 1 pixel which represents the average value of the original image z. The filter coefficients of the low pass filter of the 

20 preferred embodiment are presented in fig. 5. They correspond approximately to the samples of a two-dimensional 
gaussian distribution on a 5x5 grid. The same filter coefficients are used for the low pass filters 41 ,41\...,41 n ' at all 

scales. The same filter kernel with all coefficients multiplied by 4 is also used within the interpolators 42,42' 42"'. The 

factor 4 compensates for the insertion of zero pixel columns and rows as explained above. 

A decomposition of the original image z consisting of a background g and a residual component r 

25 

z(iJ) = 9V.j)+ r(i.j) 

is depicted in fig. 4. Each of g and r is obtained by reconstruction using the appropriate subset of pyramid levels, as 
explained next. 

30 The background image g is obtained by reconstruction using the lowest resolution levels. The coarsest resolution 
image g L is first interpolated by interpolator 44 to twice its original size and the interpolated image is next pixelwise 
added to the image of the next coarsest level r LA using adder 45. The resulting image is interpolated and added to the 

next finer detail image. If this process is iterated L - M times using the unmodified resolution levels r LA r M , then a low 

resolution background image g M will result. This image is blown up to the size of the original image z by means of the 

35 interpolators 46,46* ,46". to yield background image g. If M = L , and thedetail images //__-,,. ..,r M are left unmodified, the 
original image z will result by this reconstruction process. 

The residual or error image is obtained by reconstruction using the highest resolution levels. Starting from resolution 
level M , the detail image r MA is interpolated by interpolator 47 to twice its original size and the interpolated image is 
next pixelwise added to the image of the next finer detail level r M _ 2 in adder 48. The resulting image is interpolated and 

40 added to the next finer detail image. If this process is iterated M times using the unmodified resolution levels r M ^ r 0 , 

then a full resolution residual or error image r is obtained. 

Obviously, r and g need not be computed simultaneously: either g is computed and r can be obtained in a straight- 
forward manner as r = z - g , or r is computed and g can be obtained as g = z - r by image subtraction. 

45 4.2. Modeling and detec tion of the defects 

To the purpose of detecting defects on a recording member, the recording member is exposed to a flat field expo- 
sure, which ideally yield a constant signal amplitude referred to as the background image. The defects can be modelled 
as artefact-like shaped signals superimposed on the background variation immersed in noise. The point defects are 
so peaks of small support in both directions. Their amplitude exceeds the peak-to-peak range of the noise. Line defects 
coming from scratches and traces have one of their dimensions substantially larger than the other. Their cross-profile 
exhibits a peak-like form. Area defects introduce global deviation from a constant signal and represent variations with 
wider extent and therefore need global correction. 

Since point and line defects have at least one subtantially smaller dimension, and hence have low area, the background 
55 image is only slightly affected by the presence of the defects. This property enables the computed background approx- 
imation to pass through the defects, without actually being disturbed by them, and hence point and level defects pro- 
duce sufficient residual level signal. Area defects on the other hand do affect the parameters of the analytical 
approximation methods and distort the shape of the surface obtained by discrete or analytical approximation methods, 
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hence the calibration matrices derived will correct for them in the actual image, what should be achieved. 

4.2. 1. Local defect detection 

5 In view of the fact that the background variation is not constant, global thresholding is not straightforwardly applica- 

ble. The threshold should be floating as a function of background variation, or alternatively accomplished after back- 
ground removal. Irrespective of which of the former methods is used to approximate the background, background 
removal is generally accomplished by the following equation to yield the residual signal: 

10 r(ij)=z(i t j)-z(i,j) 

in which z(i, j) is a background fitted signal obtained by one of the formerly described methods. A defectless error signal 
has zero mean u. f and known noise variance a f . Peaks resulting from defects will be superimposed and can be detected 
by thresholding the absolute value \r(i, j)\ with a noise dependent threshold T 

15 

e(/.y)=1if|r(i,D|> T 

=0 elsewhere 

T = aa r 

20 

and a typically in the range 3. .5. For gaussianly distributed noise A/(^,ct), 99.9% of all gray values will lie in the 3-a range 
around the mean value ji. Any deviation larger than the noise-dependent threshold is attributed to defects. 

4.2.2. Global defect correction 

25 

The availability of either an analytical model or an empirically derived background surface enables the establish- 
ment of a two-dimensional calibration matrix to correct for global distortions. By global distortion is meant any deviation 
on the recording means with sufficiently large extent w.r.t. the image dimensions. This extent relates to the size of the 
smoothing kernel for empirical background fitting methods and to the shape of the basis functions for the analytical 
30 approximation methods. The correction is done in either an additive or a multiplicative manner. An additive calibration 
matrix C a (/, j) (see Figure 6) is constructed by computing the pixelwise offset between the smooth background image 
9di> i) and the global mean value \i R of the reference image obtained by a flat field exposure of the recording member 

C a VJ) = iL R -g R (i.j) 

35 

Alternatively, a multiplicative correction matrix C m (/, j) is computed (see Figure 6) as the ratio of mean value \i R 
and the smooth background variation g R (/, j) 

Application of the calibration matrices to an actual exposure as shown in Figure 7 will undo any global variation 
w.r.t. the mean value due to recording member deficiencies. 

45 Calibration of the actual image using the calibration matrix requires both images be registered. Registration of two 
images taken at different exposure and read-out times and representing physically different objects, requires the use of 
landmarks or fiducial markers during exposure. We have used prior art for the registration problem such as reported in 
Y. Christophe, J. Cornells, P. De Muynck, Subpixel geometric correction of digital images, Signal Processing IV, Proc. 
of EUSCIPCO 1988, Vol. 3, pp 1605-1608. Detection and location to subpixel accuracy of landmarks, enables the 

so establishment of a coordinate system for the actual and the reference image. An affine transformation matrix is com- 
puted, transforming the coordinate system of the actual image to the coordinate system of the reference or calibration 
image. Next, the calibration image and diagnostic image are brought into correspondence by applying the affine trans- 
formation to the actual image z(/, j) to yield the geometrically transformed image z'(/, j) . Pixel values at non-integer 
coordinates are obtained by bi-linear interpolation. Other interpolation schemes could be envisaged.. 

55 A global additive correction is performed by adding to the registered image z\U j) the correction matrix C a (/, j) or 
alternatively stated, by adding a pixel dependent offset C a (i t j) to the image z'(/, y), with R denoting the reference or 
calibration image: 
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c a </, J) = zV, J) + C a (/. j) = *'('. /) + (Mr - 9*(i, /))■ 



The multiplicative correction is performed by multiplying the registered image z'(', /) wrtn a P ixel dependent correc- 
tion factor C m (i, j) 



4.3. Postprocessing 

The result local defect detection by thresholding the residual or error image is a binary image, in which 1 denotes 
a defect (foreground) pixel, 0 denotes a background pixel. Spurious false alarms with small, user defined size can be 
removed by successive erosions of the foreground. The original thickness or remaining defects can be restored by an 
equal amount of dilations. The shape of the planar mask used was a 3x3 square mask. 

Other morphological processing could be envisaged such as presented in A.K. Jain, Fundamentals of Digital Image 
Processing, Prentice Hall, 1989, p. 384 to preserve specific shapes of defects such as lines. 

4.4. Overlay construction 

For the purpose of signalling critical areas to the operator, requiring his specific attention, an overlay image (see 
Figure 6) is constructed, representing possible defective zones detected in the previous stages. The overlay image rep- 
resents the defective zones either as coloured regions or as countours. 

4.4. Characterization and classification of the defects 

The availability of a bitmap image representing the defective pixels enables characterization and automatic classi- 
fication of the defects as either point-like, line-like or area like. For the purpose of shape characterization, shape fea- 
tures have been computed such as documented in (1) A.K. Jain, Fundamentals of image processing, Prentice Hall, 
1989, p. 390, Shape features; (2) W.K. Pratt, Digital Image Processing, 2nd ed., Wiley, p. 629, Shape analysis. Features 
of each of the connected components are computed. Geometric features include perimeter, area, minimum and maxi- 
mum radii, Euler number, corners on curved like defects, compactness, symmetry. Moment based features include 
center of mass, orientation, bounding rectangle, best fit ellipses and eccentricity. Features based on Fourier descriptors 
of the contour are derived from the amplitudes and phases of the Fourier coefficients. Other descriptors could be envis- 
aged as well. Classification techniques well known in the art such as statistical or neural classifiers have been used to 
classify a defective region as either a point, a line or an area defect. The severity of a potential defect is further assessed 
by examining the ratio of the amplitude of the defect signal with respect to the local background gray value. Based on 
shape and amplitude of the defects, the recording member finally is either definitely accepted for use, marginally 
accepted or rejected. 



1 . Method of correcting a radiation image stored in a recording member for recording member defects comprising : 

generating a calibration matrix for a recording member, said calibration matrix comprising for each pixel on said 
recording member a correction value that represents the deviation of the value that is effectively obtained in a 
pixel and a value that would he expected, when said recording medium is subjected to a flat field exposure and 
subsequently read-out by means of a calibrated read-out system, 
- generating a digital signal representation of the radiation image by reading the image by means of a calibrated 
read-out system, 

correcting each pixel of said radiation image by applying to a pixel a corresponding value of said calibration 
matrix. 

2. A method according to claim 1 wherein 

(i) said calibration matrix is obtained by 




Claims 



exposing said recording member to a flat field exposure of radiation, 
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reading a flat field image stored in the recording member exposed to said flat field exposure, 
calculating a global mean value n R by averaging pixel values of said flat field image, 

- generating a smooth background image g r (i,j) 

generating an additive calibration matrix C a (i,j) as u R - g R (i,j). and 

5 

(ii) each pixel of said radiation image is corrected by adding to a pixel a corresponding value of the additive cal- 
ibration matrix. 

3. A method according to claim 1 wherein 

10 

(i) said calibration matrix is obtained by 

exposing said recording member to a flat field exposure of radiation, 

reading a flat field image stored in the recording member upon said flat field exposure, 

- calculating a global mean value n R by averaging pixel values of said flat field image, 
generating a smooth background image g r (i,j) 

= generating amuttiplirativec^ 

(ii) each pixel of said radiation image is corrected by multiplying a pixel with a corresponding value of the mul- 
20 tiplicative calibration matrix. 

4. Method of localising local defects in a recording member for recording a radiation image comprising the steps of 

exposing said recording member to a flat field exposure of radiation, 
25 - reading a flat field image stored in the recording member upon said flat field exposure, 

generating a smooth background image g R (i,j) 

generating a residual image r(ij) by subtracting from a pixel value of the flat field image a corresponding value 
of the smooth background image g R (i,j). 
- thresholding the absolute value |r(i,j)| of said residual image r(i,j) with a noise dependent threshold T to gener- 
30 ate an error image e(ij) wherein 

e(i.j) = 1 if |r(i,j)| > T, 
e(i,j) = 0 elsewhere, 

35 

and wherein T = aa R , a being comprised between 3 and 5 and 

a R heing the variance of the noise distribution in said flat field image 

5. A method according to any of claims 2 to 4 wherein said smooth background image g R (i J) is obtained by applying 
40 a regression technique to pixel values of said flat field image. 

6. A method according to any of the claims 2 to 4 wherein said smooth background image g R (i,j) is obtained by apply- 
ing a spline fitting technique to pixel values of said flat field image. 

45 7. A method according to any of the claims 2 to 4 wherein said smooth background image g R (i,j) is obtained by 
smoothing pixel values of said flat field image. 

8. A method according to claims 2 or 4 wherein said digital signal representation of original image 'z* is decomposed 
into a sequence of detail images r 1t r 2 ... r^ and a rest image g L at a still lower resolution level and wherein said 

so residual image Y is formed by summing up detail images at M finest resolution levels starting from a level M - 1 up 

to level 0, each of the detail images being interpolated to the size of a detail image at a next finer resolution level 
before being summed up. 

9. A method according to claims 2 or 4 wherein said digital signal respresentation of original image *z' is decomposed 
55 into a sequence of detail images r 1( r 2 , ... r u1 and a rest image g L at a still lower resolution level and wherein said 

residual image V is formed by 

(i) summing the rest image g L and the detail images r { at L-M coarsest resolution levels, starting from resolution 
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level L up to resolution level M, each being interpolated to the size of a detail image at a next finer resolution 
being summed, to form a downzoomed background image g M , 

(ii) upzooming background image g M to form a full-size background image 'g* by performing M successive inter- 
polations on g M , 

5 - (iii) subtracting the upzoomed background image 'g* from the digital signal representation of original image *z\ 

10. A method according to claim 1 wherein said calibration matrix is generated in advance, stored in memory and 
retrieved upon identification of a recording member. 

io 11. A method according to claim 4 wherein said error image e(i,j) is applied to a hard copy recorder to generate an 
overlay image. 

12. A method according to claim 4 wherein said error image e(i,j) is subjected to characterisation and classification to 
generate a report on shape and amplitude of defects. 

15 

13. A method according to claim 1 wherein (i) said recording member is a photostimulable phosphor screen and 
wherein (ii) a radiation image is read-out by scanning said screen by means of stimulating radiation, detecting light 
emitted upon stimulation, converting detected light into an electric signal and digitizing said signal. 
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